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Abstract 

We study the statistical mechanics of classical Coulomb systems in a 
low coupling regime (Debye-Hiickel regime) in a confined geometry with 
Dirichlet boundary conditions. We use a method recently developed by the 
authors which relates the grand partition function of a Coulomb system 
in a confined geometry with a certain regularization of the determinant of 
the Laplacian on that geometry with Dirichlet boundary conditions. We 
study several examples of confining geometry in two and three dimensions 
and semi-confined geometries where the system is confined only in one or 
two directions of the space. We also generalize the method to study sys- 
tems confined in arbitrary geometries with smooth boundary. We find 
a relation between the expansion for small argument of the heat kernel 
of the Laplacian and the large-size expansion of the grand potential of 
the Coulomb system. This allow us to find the finite-size expansion of the 
grand potential of the system in general. We recover known results for the 
bulk grand potential (in two and three dimensions) and the surface ten- 
sion (for two dimensional systems). We find the surface tension for three 
dimensional systems. For two dimensional systems our general calculation 
of the finite-size expansion gives a proof of the existence a universal loga- 
rithmic finite-size correction predicted some time ago, at least in the low 
coupling regime. For three dimensional systems we obtain a prediction 
for the perimeter correction to the grand potential of a confined system. 

I Introduction 

The study of confined classical Coulomb systems have attracted attention for 
some time in particular because in some cases they exhibit universal proper- 
ties J] |21 El- This universal behavior is present in some correlations functions 4 
and also on the thermodynamic quantities of the Coulomb systems. In particu- 
lar, the grand potential and thus the free energy, exhibit finite-size corrections 
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which depend on the confining domain 1 . Two particular cases have attracted 
attention: fully confined systems and semi-confined systems. On this respect we 
shall speak of semi-infinite or semi-confined systems to refer to systems which 
are confined only in a certain spatial direction (for example systems confined in 
a slab), to distinguish them from the totally confined systems. In both cases, 
exactly solvable models in two dimensions have allowed the explicit calculation 
of the finite-size corrections in the free energy for a given value of the coupling 
constant. For systems confined in a metallic slab of width W, in d dimensions, 
the free energy and the grand potential per unit area (times the inverse reduced 
temperature) exhibit an algebraic universal correction C(d)/W d ~ 1 with 

OM-^ (...) 

where T(z) and £(z) are the Gamma function and the Ricmann zeta function 
respectively. This has been shown £Q to hold for any general Coulomb system 
provided that the system is in a conducting phase and it has good screening 
properties. It has also been checked in several solvable models. The correction 
is universal in the sense that it does not depend on the details of the microscopic 
constitution of the system. 

For two-dimensional fully confined Coulomb systems there are also univer- 
sal finite-size corrections which are similar to those of two dimensional critical 
systems. This lead us to another interesting feature of conducting classical 
Coulomb systems, which is its manifest similarity with critical systems 013]. 
Although the particle and charge correlation functions of the Coulomb system 
are short-ranged because of the screening, it has been shown that the correla- 
tions of the electric field and the electric potential are long ranged @JE]- In this 
sense they can be considered as critical systems and they share properties of 
statistical models at criticality. For example, in two dimensions conformal field 
theory, which has proved to describe and classify correctly critical systems, pre- 
dicts the existence of universal corrections in the free energy for critical systems 
due to their finite-size jSJ [7| |S] • Explicitly, for any two-dimensional statistical 
system in its critical point confined in a domain of characteristic size R with 
smooth boundary, the free energy F has a large- R expansion of the form [SI {7\ 

0F = AR 2 + BR - — In R + ■ ■ • (1.2) 
6 

where f3 — X/iksT) with T the absolute temperature and ks the Boltzmann 
constant. The first two terms AR 2 and BR represent respectively the bulk free 
energy and the "surface" (perimeter in two dimensions) contribution to the free 
energy. In general, the coefficients A and B are non-universal (they depend on 
the microscopic detail of the model under consideration) but the dimensionlcss 
coefficient of In R is highly universal depending only on the Euler characteristic 
of the manifold \ = 2 — 2h — b, where h is the number of handles and b is the 
number of boundaries, and on c the central charge of the model. For Coulomb 
systems the existence of a similar expansion, which reads 

(3F = AR 2 +BR+-\nR+--- (1.3) 
6 
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has been shown to hold in several exactly solvable models at a fixed value of the 
coulombic coupling constant ^ [31 El HOI ITT] and in some particular geometries 
for any value of the coupling ^] ED ^] E] • 

In a previous paper |16) we considered two-dimensional Coulomb systems in 
a low coupling regime, the Debye-Huckel regime. We computed the grand po- 
tential for systems confined in two simple geometries, the disk and the annulus 
with ideal conductor boundaries, and confirmed the validity of the finite-size 
expansion i|1.3[) in those cases. We showed that the grand canonical partition 
function for a classical Coulomb system in the Debye-Hiickcl regime, confined 
with ideal grounded conductor boundaries, can be expressed as an infinite prod- 
uct of functions of the eigenvalues of the Laplace operator satisfying Dirichlct 
boundary conditions. The explicit form of this spectrum and the corresponding 
infinite products, depend on the shape of the confining domain, and must be 
calculated for each particular geometry. By a careful calculation of these infi- 
nite products we obtained the explicit form of the grand potential for Coulomb 
systems confined in a disk and in an annulus. When these systems are large 
we computed the finite-size expansion of the grand potential and we found the 
universal correction predicted by Eq. I|1.3|> . The first purpose of the present 
paper is to apply this method to other particular cases of confining geometry 
including semi-confined systems and also to systems in three dimensions, for 
which conformal field theory predictions do not apply. 

The second purpose this paper has to do with the fact that, from a more 
general point of view, it is possible to define a spectral function for the Lapla- 
cian, the heat kernel, that turns out to have an asymptotic behavior for small 
argument which is independent of the explicit form of the eigenvalues |17l I18| . 
Making use of these results, it is possible to show that the spectrum of the 
Laplace operator calculated on a given manifold endowed with a metric, contain 
geometrical information about the manifold itself. In this paper we use those 
ideas to obtain the large-size expansion of the grand potential for Coulomb sys- 
tems confined in arbitrary geometries. Our results for the particular cases agree 
with the predictions of this general formalism. 

This paper is organized as follows. In section^ we summarize a few results 
of our previous paper \nj\ concerning the calculation of the grand potential for 
Coulomb systems in the Debye-Huckel regime in given confining geometries. In 
particular, we briefly describe how the grand potential can be obtained in terms 
of an infinite product of functions of the eigenvalues of the Laplace operator. 
In sections IlIII and llVl we apply the general method from Ref. ^Hj reviewed in 
section m In section ITTTI we apply the method to some particular examples of 
confined and semi-confined systems in two and three dimensions. In section llVl 
we consider the general case of fully confined systems in an arbitrary geome- 
try. We relate the grand potential of the system to the zeta regularization of 
the determinant of the Laplacian. By using the known results ^| EH for the 
asymptotic expansion of the heat kernel we find in general the finite-size expan- 
sion of the grand potential and, for two-dimensional systems, we confirm the 
existence of the predicted universal finite-size expansion. At the end of that 
section we present an illustration of this latter method by considering the case 
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of a Coulomb system confined in a large square, and we recover a finite-size 
correction predicted by conformal field theory. Sections II I II and IfVI arc mostly 
independent and the reader not interested in the examples of section 11111 can 
proceed directly to the general treatment exposed in section llVl In section Ivl 
we present a summary and gather some conclusions. 



II Summary of Previous Results 

Let us start by describing the model under consideration. Our system is a multi- 
component Coulomb gas living in d dimensions and composed of s species of 
charged particles a = 1, s each of which have N a particles of charge q a . The 
system is confined in a domain of volume V with ideal conductor boundaries. 
We shall describe describe this systems using classical (i.e. non-quantum) sta- 
tistical mechanics in the grand canonical ensemble. The average densities of the 
particles n a are therefore controlled by the fugacities £ a . We shall impose the 
pseudo-neutrality condition 

J>aC Q =0 (2.1) 

a 

which implies that at the mean field level the system is neutral and there is no 
potential difference between the system and the boundaries. In the appendix 
B of Ref. |1K | we explain what happens in the more general case when the 
condition l|2.1|l is not satisfied. 

The interaction potential between two unit charges located at r and r' is 
given by the Coulomb potential v(r, r') which is the solution of Poisson equation 

Au(r,r') = -s d S(r-r') (2.2) 

satisfying Dirichlet boundary conditions and where Sd = 27r d / 2 /r(<i/2), that is 
in two dimensions S2 = 2-7T and for three dimensional systems S3 = An. For 
non-confined systems the Coulomb potential reads 

' if d 



v°(r,r') = I Ir-r'l (2.3) 



- In t if d = 2 

where L is an arbitrary length scale which fixes the zero of the Coulomb potential 
in two dimensions. For the confined system under consideration the explicit form 
of the Coulomb potential should be modified in order to satisfy the Dirichlet 
boundary conditions. 

As explained in Ref. ^Hl ( see a ls° EHX the potential energy of the system 
can be written as 



= I E E 9«97^( r «^ r 7j) - 1 1! E 9 «"'( rtt ." r ^') ( 2 - 4 ) 
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In the first line the prime in the first summation means that the case when 
a = 7 and i = j must be omitted. The first term is the usual inter-particle 
energy between pairs. The second term is the Coulomb energy of a particle 
and the polarization surface charge density that the particle has induced in the 
boundaries of the system. 

In Ref. |16j we performed the sine-Gordon transformation |2(J| on the grand 
canonical partition function 5 of the system. Then we expanded the action 
around the mean field solution to the quadratic order in the field. This is valid 
in the Debye-Hiickcl regime. Then the remaining functional integral can be 
performed easily since it is a Gaussian. The result is a certain determinant 
involving the Laplacian which we put in the form 

where n" 1 — (J2 a SdCa0Qa) equals the Debye length in this regime, A m 
denotes the Laplacian eigenvalues satisfying the Dirichlet boundary conditions 
and A° = — K 2 , K € M. d , refers to the (continuum) eigenvalues of the Laplacian 
in the non-confined case. These come from the "subtraction" of the self-energy 
term w (r,r) in Eq. (|2.4|) . 

Each infinite product in (|2.5[) diverges separately. Indeed they are ultravio- 
let divergent for large values of |A m | and |A°|. However, when they are putted 
together as in l|2.5|) . the divergences cancel each other (at least for the bulk 
properties of the system). In three dimensions we can find immediately a well 
defined expression for the grand potential from il = — fc^TlnS. In two di- 
mensions, the situation is a bit more involved since certain infrared divergence 
appear in the second product and it must be regularized by introducing a lower 
cutoff. In Ref. ^Sj we explained how to deal with this case and we found the 
value of this cutoff explicitly in terms of the constant L which fixes the zero 
of the Coulomb potential, which needs to be supposed to be large. This cutoff 
was found to be given by fc min = 2e~ c /L where C is the Euler constant. From 
Ref. JS| j we recall that for a non-confined system Eq. (|2.5|l gives for the bulk 
grand potential 



, kL ~ 1 
-In C + - 

2 2 



Co (2D) (2.6) 

a 

J2(a (m (2.7) 



in two and three dimensions respectively. These expressions agree with results 
by the usual formulation of the Debye-Hiickel theory G3] • 
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Ill Solved Examples 



III.l Systems in Two Dimensions 

In Ref. IBj using Eq. (|2.5II we computed the grand potential of a two-dimensional 
Coulomb system confined in a disk and in an annulus, and we confirmed that 
its finite-size expansion is of the form 

pn= f3n b + (3-f+^\n( K R) + 0(l) (3.1) 
6 

with Ofc, given by Eq. (|2.6(l . The surface (perimeter) tension is 

7 = -fc B T K /8. (3.2) 

We notice the existence of the universal finite-size correction (x/6) In R with the 
Euler characteristic x ~ 1 for the disk and x = for the annulus. 

In this section we will consider an additional example of confining geometry 
in two dimensions. 



III. 1.1 Space Between Two Infinite Lines: The Slab in Two Dimen- 
sions 

The method outlined in section [D] can be used to study semi-confined systems. 
In this subsection we consider the case of such a system in two dimensions. 
The geometry consists of two infinite parallel lines spaced by a distance W 
and the Coulomb systems is confined in between these two lines. We assume 
Dirichlet boundary conditions for the electric potential. Let us assume that 
the lines are in the direction of the y-axis and the z-axis is perpendicular to 
the lines. If we write the Laplacian eigenvalues as A = —k\ — ky, these take 
discrete values only in the k x direction. The eigenfunctions can be written as 
fy(x,y) oc e l ( kyy } sin(k x x), satisfying the boundary conditions ^(0,y) — and 
fy(W, y) = 0, which imply k x = nir/W with n a positive, non-zero, integer. In 
the direction of the y-axis there is no confinement therefore k y € R. Then, the 
eigenvalues of the Laplace operator are given by \ n ,k y — ~ (mr/W) 2 — ky, for 
n = 1, 2, . . . , and k y S K. Introducing the explicit form of the eigenvalues in 
(|2.5|) we have the grand potential expressed as 

-i 7 poo 00 / 2 \ 1 2 

- w L H (i + <^j) *. - 5 e k - ? <"> 

where I is the dimension of the system in the y-direction. The second term 
in H3.3(l involve the spectrum for a non-confined system, A^ -* = — K 2 with 
Kel 2 . It can be written as \ J2 k ^ = Jj*™* ^f, where V — IW is the 

"volume" (area) of the system between a portion of length I of the confining lines. 
The lower limit for this integral is k m i n — 2e~ c / L as mentioned in section HTl 
and explained in Ref. [lfi| . Also as explained earlier this integral is ultraviolet 
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divergent therefore we introduced an ultraviolet cutoff K max . In the first term 
of Eq. (|3.3[1 the infinite product converges to a known expression [21] giving 



(w^Jk 2 + kfj 



l k y Sillh [Wjtf + k*) mK 2 f K^ dR 



1 ■'" /- 2 1 fc^ sinh(fc !/ ^) 4?r •/*»*. ^ „ 

(3-4) 

where we also introduced an ultraviolet cutoff for the first integral fc max . Both 
should be proportional and their exact relation can be 
found by requiring that in the limit W — > oo we recover the known bulk value 
of the grand potential (|2.6|) . Performing some of the integrals in (|3.4|) we find 
that the grand potential per unit length w = fi/Z is given by 



K 2 W^ 2fc max K 2 W 

P w = — ; — l n ~~t? 1 — ; — 

47T A ma .x 47T 



I 

t: — hi tt - 

2 2e" 



c 



Wj^Cc (3.5) 
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where all terms that vanish when fc max — * oo have been omitted. Therefore to 
recover the known value (|2.6|) of the bulk grand potential in the limit W — > oo 
the ultraviolet cutoffs should be related by -Km ax = 2fc max . Using these cutoffs 
we finally find the grand potential and its finite-size expansion 

fr> = /M + 2/3 7+ ^ 7 + ^ 00 ln(l- e -^V^) ^ (3.6a) 

= fa b + 2(3 7 +^ + 0(e- 2 « w ) (3.6b) 

where = fib /I with the bulk grand potential fif, given by Eq. (|2.6|l . The 
surface tension 7 is given by Eq. I|3.2|l . which is the same surface tension that 
we found in the case of the disk and the annulus in our previous work |16j as 
expected. Finally, we also found the universal finite-size correction for the case 
of the slab in two dimensions, which turns out to be n/(24W) = £(2)/(47rW) in 
accordance with the general prediction from Ref. [Q: T(d/2)((d)/(2 d Tr d / 2 W d ~ 1 ) 
for d = 2 as expected. 

III. 2 Systems in Three Dimensions 

In this section we consider some examples of three-dimensional Coulomb systems 
first confined in a slab geometry then inside a ball and inside a spherical thick 
shell. 

III. 2.1 Space Between Two Infinite Planes: The Slab in Three Di- 
mensions 

We begin the study of particular examples of three-dimensional systems by 
considering a system confined in the space between two infinite parallel planes, 
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separated by a distance W . Taking the ^-coordinate along the direction normal 
to the planes we find that the eigenfunctions are ^(r) oc e l ( k ±- r ±) sin(k x x), where 
kj^ ■ rj_ = yk y + zk z , and satisfying the boundary conditions ^(0,y, z) = = 
&(W, y, z). Thus, the eigenvalues are given by k x = nir/W, with n = 1,2,..., 
and k y € M and k z e R. Using l|2.5[l and the explicit form of the eigenvalues we 
have 



0Q 



1 A 



2 (2tt) 



mn 

n=l 



iw) + k i, 



dkj_ 



2^X°, 



VC a (3.7) 



where A represents the area of the planes. The second term coming from the 
subtraction of the self-energy is now 



1 ^ k 2 
2^ A° 



1 K 2 V 

2 (2tt) 3 



dn 



K„ 



k 2 dk 

k 2 



K V 

(2tt): 



(3.8) 



where V = AW is the volume of the system. As in the two-dimensional example 
we introduced an ultraviolet cutoff K max . 

Similarly to the two-dimensional slab, the infinite product in the first term 
converges to [21] 



n 



sinh [W^/k 2 + k 



( H2L\ 
V W ) 



V^ 2 



sinh{k±_W) 



(3.9) 



with k± = | k_i_ | . The remaining integral over kj_ is ultraviolet divergent and 
must be cutoff to a maximum value A: max 
have 



for k i . In the limit fc n 



oo we 



~A 



K 

8^ 

~8n 



2Kr, 



In- 



C(3) 
16irW 2 



127T 



i r 00 

— / fclnfl 



D -2WVk' 2 +K^ 



)dk 



(3.10) 

(3.11) 
(3.12) 



From this expression we see that the ultraviolet cutoffs should be related by 
fcmax = 2K max /Tr in order to recover the known value of the bulk grand poten- 
tial (|2.7|l in the limit W — > oo. Then the grand potential is finally given by 



v ; „ H / uln(l 

I671W 2 4tt J 1 y 



-2ukW 



A 



2h- 



C(3) 
I&ttW 2 



0{e 



-2kW\ 



)du (3.13a) 
(3.13b) 
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with the bulk grand potential Qb given by Eq. i|2.7[) and the surface tension 7 
given by 

1 



1D7T 



In — 



H 2 



w ln I— ( 3 - 14 ) 

00 lb7T /Cmax 



Note that when we take the limit fc max — > 00 the surface tension diverges with 
the cutoff as — [k 2 /(167t)] lnfc max . This divergence in the surface tension can be 
understood if we note that the particles tend to move to the frontiers because 
of the ideal conductor character of the boundaries. This is easy to see from a 
physical argument: the ideal conductor boundaries condition is equivalent to 
introduce an image charge of opposite sign at the other side of the boundary for 
each particle in the system. Particles near the boundary "feel" an attraction to 
the boundary due to their proximity with their corresponding images. Near a 
boundary, the density of the species a at a distance X from the boundary will 
behave in this low coupling approximation as the linearized Boltzmann factor of 
the particle-image interaction 1 + f3q 2 t /(4X). At large distances this interaction 
is screened but at short distances it remains non-integrable. Since the surface 
tension can be obtained as an integral of the density profile [25], this surface 
tension will be infinite. Imposing a short-distance cutoff D for the minimum ap- 
proach of the particles to the wall, will give a the surface tension which diverges 
as In IX Our ultraviolet cutoff k max is proportional to 1/D. For details see 
Ref. (T^|. Notice that, on the other hand, for a two-dimensional system the sur- 
face tension does not diverge with the cutoff (see Ref. ^5] and section UlL 1. 1|) . 
In two dimensions the particle-image interaction is [<z 2 /2] ln(2X/L) and this 
expression is integrable at short-distances. This explains why the surface ten- 
sion is finite and cutoff independent for two-dimensional systems although the 
particles are strongly attracted to the boundaries contrary to the situation in 
three dimensions where the surface tension diverges with the cutoff. 

Returning to Eq. i|3.13bjl we found a finite-size correction depending on W~ 2 . 
This agree with the universal finite-size correction for a slab in <i-dimensions, 
Eq. (|I.1|) . for d = 3 predicted in Ref. [T]. 



III. 2. 2 Coulomb System Inside a Ball 

We continue the study of finite-size Coulomb systems calculating the grand 
potential for a three-dimensional Coulomb system confined inside a spherical 
domain. The eigenvalue problem for the Laplace operator in this case is easily 
solved. The eigenf unctions are ^(r,9,ip) = \J n/(2rX 1 / 2 )Ii + i/ 2 (VXr)Yi m (9, tp) 
where h+1/2 are the modified Bessel functions of half integer order and Yi m are 
the spherical harmonics. The eigenvalues A are determined from the Dirichlct 
boundary condition ^(R, 0, (f) = where R is the radius of the sphere. Thus, the 
eigenvalues are the roots of the equation Ii + i/2(VAi?) = 0. Let us call z^+1/2, n 
the zeros of Ii+\/ 2 - Then the eigenvalues are given by Afc = v 2 +l / 2 f° r 

I = 0, 1, 2 . . . and n = 1, 2, Also for each value of I and n, the corresponding 

eigenvalue is degenerated 21 + 1 times. Then, the expression for the grand 
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potential obtained from l|2.5[l takes the form 



In 



21+1 



nn> 

1=0 \n=l \ 



1+1/2,' 



^Ew-E^ (3.15) 



where the indexes n and I + 1/2 denote the root and the order of the modified 
Bessel function I1+1/2 respectively and z — kR. The infinite product over the 
index n can be performed exactly |24j 



n 



1 - 



Z+l/2,n / 



= r i+ - 



1+1/2 



h+i/2{z). 



(3.16) 



The remaining summation over the index / diverges and we must regularize it 
by introducing an upper cutoff N on I. 

The second term coming from the subtraction of the self energy is 



1 K 2 V 

2 (27t) 3 



(19. 



k 2 dk 



k 2 



— R- 3 K 

— o il J1 max 
6TT 



(3.17) 



where V = \^R? is the volume of the system. As in the previous examples we 
introduced an ultraviolet cutoff i^ max which must be proportional to N in order 
to cancel the divergences. The exact relation between K mayi and N is found by 
the requirement that the bulk value of the grand potential (|2.7|l is recovered in 
the limit R — > 00. 

To find the finite-size expansion of the grand potential we make use of the 
Debye uniform asymptotic expansion [2E] for the Bessel functions, valid for large 
values of the argument, 



lnl u (z) = -iln(27r)-iln(z 2 + j/ 2 )+?7(^,z) + - 3u 



24^ 



1 



z 2 + v 1 



with 



rj{v,z) = (z 2 + z/ 2 ) 



2\l/2 



v sinh 



[z 2 + u 2 ) 1/2 



(3.18) 



(3.19) 



together with the Euler-McLaurin summation formula to transform the summa- 
tion into an integral: f(l) = f" f{l)dl+\ [/(0) + f(N)]+± [f(N) - f(0)} + 

• • • , and the Stirling asymptotic expansion for the Gamma function. Then in 
the limit ./V — > 00 and z — > 00 we find 



1 + 2 In 

N 



Rk\ k 2 R 2 kR 



8 



^ + O(R )+o(N°) 



where 



/5^buii 



K 

127 



3k 2 N 
16^ ~R 



4tt 2 



V 



(3.20) 



(3.21) 
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This value of the bulk grand potential should be equal to the one given by 
Eq. (|2.7|l therefore K max = (3w/4)N/R. The ultraviolet cutoff K max is indeed 
proportional to the cutoff N and with this relation between the cutoffs the bulk 
divergences of the first and second terms in the r.h.s. of l|3.15(l cancel each other. 
Reporting the result for if max in 13.21[1 we find 

(3.22) 

The first term in the r.h.s. of Eq. I|3.22|l is the bulk grand potential in three 
dimensions. In the second term, the expression in parenthesis gives the surface 
tension for the system and the third term is a (non-universal) perimeter contri- 
bution. Notice again that the surface tension diverges as — [k 2 /(167t)] \nK max 
with the cutoff. 



0Q = 



12tt 



III.2.3 Spherical Shell 



We consider now the case of a three-dimensional Coulomb system confined inside 
a spherical shell. Let a and b be the inner and outer radii of the shell respectively. 
As in the previous example, the eigenfunctions ^S(r, 9, tp) are separable in a radial 
and an angular part. The eigenfunctions of the Laplacian for this geometry are 



*(r,0,<P) = A^=I l+1/2 (^Xr) + B^=K l+1/2 



(VXr 



Y lm (0,<p). The 

eigenvalues are determined by the boundary conditions W(a, 9, ip) — \&(&, 9, if) = 
0, that is AI l+1/2 (^/Xa) + BKi +1/2 (y/Xa) = = AI l+1/3 (V\b) + BK l+1/2 (V\b). 
These two equations can be considered as a linear system of equations for the 
coefficients A and B. It has a non-trivial solution if and only if 



4+1/2 ( VXa)K l+1/2 ( VXb) - Ii +1/2 (V\a)K l+1/2 (VXb) = . 



(3.23) 



The roots of this equation gives the eigenvalues. Let i?; iTS be the n-th root of 

Ii+i/ 2 {za)K l+1/2 (zb) - Ii+i/ 2 {za)K l+1/2 (zb) = for I = 0, 1, 2 Then we 

have Afe = $f n . Each eigenvalue is (21 + l)-degenerated. Then the expression 
for the grand potential takes the form 



i 



■In 



21+1 



nn 

1=0 n=l 



,q2 



(3.24) 



The infinite product over the index n can be performed explicitly by using a 
method explained in Refs. |27M10llTS] . Let us consider the entire function 



/,(*) = (2/ + 1) 



Ii+i/2{za)K l+1/2 (zb) - Ii+i/ 2 (za)K l+1/2 (zb) 



(!) 



(!) 



(3.25) 



which has the following properties: fi(z) = fi(—z), fi(0) — 1, //(0) = and 
its zeros are "ff^n- Therefore it admits an expansion as a Weierstrass infinite 
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product 

M*) = nf 1 -^;) • ( 3 - 26 ) 

Then the product we need to compute is simply /;(/t). The grand potential then 
takes the form 

1 N 2 — , 

0Q = - £(21 + 1) In /,(«) - |- (& 3 - a 3 ) K max - £ V( a (3.27) 

Z=0 a 

where we introduced the now familiar cutoffs N for the sum and K m&x for the 
integral. As in all the previous examples these cutoffs are proportional and their 
exact relation is found by requiring that in the limit of a large system we recover 
the known bulk value (|2.7() of the grand potential. 

For the calculation of the finite-size expansion of the grand potential, notice 
first that the contribution of Ki{nb)Ii{Ka) is of order e~( b ~ a ' K , that is to say 
exponentially smaller than the contribution of the term Ii(Kb)Ki{na) and as a 
consequence we can ignore it. To find the value of the remaining summation 
over the index /, we use the approximations of K l+1 / 2 (na) and ^+1/2 for 
large values of the argument, which can be obtained from the Debye approxima- 
tion for the modified Bessel functions Then we apply the Euler-McLaurin 
summation formula and take the limits N — * 00, b — > 00, a — > 00 and b— a — > 00 
with a/b < 1 fixed. We find that the volume bulk part of the grand potential is 
given by 

„ Nn 3 o Nk 3 , Ik K 2 fc max 47r v—v \ 1 1 ^\ 

^ bulk = _*-_«»- ^ + ^ + y £ C. J (* 3 - a 3 ) (3.28) 

Equating this to the grand potential of the bulk system 1|2.7|) we find 

„ _3irNl-(a 2 /b 2 ) 



4 b l-{a 3 /b 3 ) 

Using this relation we finally obtain 



(3.29) 



/3Q = [3n b + Surface + \ (b - a) K + O(o°) + 0(b°) (3.30) 



3 

where fib is given by Eq. (jZH) and 



sur 



face 



i(21nf -3)aV + i 21n^ + 1 W (3.31) 



Some additional comments are in order. The divergence in the surface tension, 
familiar to us at this point, is also present in this case. In the limit K max — > oo we 
have /?n su rface —* ~47r (b 2 + a 2 ) [k 2 /167t] ln(K max / k) which allows us to define 
a surface tension 7 similar to the one of the previous other three-dimensional 
cases 

(3 7 =-J?-ln^ (3.32) 

l07T k 
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We find again a perimeter contribution to the grand which is not universal 
(it depends on the Debye length This time it is given by ksTn^b — a)/3. 

This clearly suggests that this perimeter correction is related to the curvature of 
the boundary and probably to the curvature of the space itself. This is indeed 
the case as we will show for any general geometry in the following section. 

IV Grand Potential for Arbitrary Confining Ge- 
ometries 

Up to now we have been capable to find the explicit form of the grand potential 
for systems in specified confining geometries. Our calculations of last section al- 
ways involve the resolution of the Laplacian eigenvalue problem for each specific 
geometry. From a more general point of view it is possible to define functions 
of the spectrum of the Laplacian that admit asymptotic expansions, that turn 
out to have some properties that are independent of the explicit form of the 
eigenvalues. Also in some cases these functions are related to some invariants 
of the confining manifold, for example in two dimensions to the Euler charac- 
teristic of the manifold. In this section we make use of those ideas to find the 
finite-size expansion of the grand potential in the case when the spectrum of 
the Laplacian in the confining geometry is not known explicitly. 

IV. 1 Spectral Functions of the Laplacian 

In this section we review some spectral functions that will be useful for our 
analysis. Let Ai g be a Riemannian manifold endowed with a certain metric 
g and with boundary dA4 g and A the Laplace operator defined on A4 g . The 
spectrum of M g is the set {0 ^ A ^ Ai ^ • • • J, —00} of eigenvalues of A, that 
satisfy A^> = \*ff, where \& are the eigenfunctions of A. In order to determine 
the spectrum, these functions must satisfy certain boundary conditions on dA4 g 
which we impose to be of the Dirichlet type, that is \& = on dM. g . The first 
spectral function we are interested in is the heat kernel defined as 

00 

e(t) = XV A *, (4.i) 

which is convergent for die(t) > 0. It is known that the heat kernel admits an 
asymptotic expansion for t — > of the form 

00 

e(t)^J2 c iJ in - ( 4 - 2 ) 

n=0 

Here {i n } is a certain increasing sequence of real numbers and < 0. The 
exponent iq is particularly important because Ci o i io is the divergent leading 
term in the series. According to the famous Weyl estimate |2U for the Laplacian 
io = —d/2 where d is the dimension of the manifold and (Air) d / 2 Ci is the volume 
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of the manifold. Following Ref. we define the order of the sequence of the 
eigenvalues as [i = — io = d/2. For a manifold M. g with a smooth boundary, 
some of the terms in the small-i expansion of the heat kernel have been found 
by Kac and others [T71ITS| 

(4 7 ri)' i/2 e(t) = V - + -(2C + D) + o(t 3 / 2 ) (4.3) 

4 6 

where 

V = the volume of M g (4.4) 
B = the surface area of dM. g (4-5) 

C = the curvatura integra = / K (4-6) 

JMg 

D = the integrated mean curvature = / J (4-7) 

JdMg 

where K is the scalar curvature at a point inside the domain M g and J is the 
mean curvature at a point in the boundary dA4 g . If we choose a metric g in 
which locally the first coordinate is perpendicular to the boundary and outward 
pointing to it then the mean curvature J can be computed as |XH| 1 

J = di[ 5 11 det 5 ]Vgi7/det.g (4.8) 

Notice that in two dimensions the well-known Gauss-Bonnet theorem [20] states 
that 2C+D = 4ttx where \ is the Euler characteristic of the manifold. Therefore 
in two dimensions the heat kernel expansion reads 

The second spectral function we are interested in is the Fredholm determi- 
nant defined as Jlfc^Lo (l ~ Tl) wn i cn ^ s precisely the infinite product involved 
in the calculation of the grand potential of the Coulomb system from Eq. 12. 5|) . 
Unfortunately this infinite product only converges for sequences of order fi < 1, 
and therefore it diverges for the cases we are interested in, when [i = 1 (two 
dimensions) and /i = 3/2 (three dimensions). For the cases /i > 1 a Weierstrass 
canonical regularization of the Fredholm determinant reads |29| 



n l 1 -ir ex p 



a a 2 oM 



^ 2A 2 fc - MA M 



(4.10) 



which is valid for fj, > 1, where [[/] is the integer part of /i. The exponential 
term in H4.10|) is introduced in order to make the infinite product convergent 

1 Notice that we use here the convention of outward pointing normal vectors to the bound- 
ary. This is the opposite convention as the one in Ref. 1181 : our J is minus the J of Ref. 1181 . 
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when the order of the sequence is larger than one. We are interested in two- 
and three-dimensional manifolds when fi equals 1 or 3/2 respectively. In both 
cases expression (|4.1U|I reduces to 

OO , v 

r(a) = n( 1 -^)e a/Afc . (4.11) 



fe=0 



Although, strictly speaking the product Ij4.1t)|l is only defined for ^ > 1 we will 
still use it in the two dimensional case when /x = 1. We will see that using the 
regularization (|4.11|) for the ultraviolet divergence of the Fredholm determinant 
in the case [i = 1 introduces some infrared divergences. However, as we will 
show in detail later, these infrared divergences can be dealt in a similar way as 
it was done for the two dimensional examples of Ref . ^Sj ■ 

We finally introduce the generalized zeta function defined as 

oo 

Z(s,a) = ]T(a-A fc r s (4.12a) 

k=0 



— / efile^'r'A. (4.12b) 
(s) Jo 



r( s ) 

The first expression for the generalized zeta function defined as a series is con- 
vergent for any s such that 3?e(s) > fi and for any a such that a ^ \q. The 
second expression l|4.12bjl where the zeta function is expressed as a Mellin trans- 
form of the heat kernel actually allows an analytic continuation of Z(s 7 a) for 
5Re(s) < /i if the heat kernel admits a full asymptotic expansion for t — > 
of the form l|4.2|) as explained in Ref. Note that Z(s,0) = Z(s), where 

Z{s) = X)fe=o ( — ^k) is the zeta function of the sequence {A&}. The analytic 
continuation of the zeta function has the following properties |29) which we will 
need shortly. Z(s) is meromorphic in the whole complex s plane and has poles 
at s = —i n with residue 

ResZ(- ln ) = -^-. (4.13) 

In particular from Eq. (|4.3(l we deduce that the first pole is encountered at 
s = —i = (j, = d/2 and has residue V/[r(d/2)(47r) d / 2 ]. Notice that this residue 
is independent of the shape of the manifold: it only depends on its total volume 
V. Also the negative or zero integers s = —n are regular points of Z(s) and we 
have 

Z(-n) = (-l) n n!c_ n . (4.14) 

The generalized zeta function provides another regularization for the Fredholm 
determinant of the Laplacian known as the zeta regularization. Indeed dif- 
ferentiating H4.12a|) under the sum with respect to the variable s and putting 
s = afterwards (a procedure which is not legal since the expression Ij4.12a|l is 
convergent only for 5Re(s) > [i) formally yields 



dZ(s,0) 



ds 



dZ(s,a) 



s=0 93 



s=0 



k=0 
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Strictly speaking this Eq. (|4.15|l is incorrect and it should be only understood 
as a formal relation to justify the word "determinant" in the name zeta regular- 
ization of the determinant of the Laplacian since the infinite product involved 
in the r.h.s. is divergent. Notice however the l.h.s. of Eq. I|4.15|l is well defined 
once the analytic continuation of Z(s, a) is done with the aid of relation (|4.12b|) . 

The zeta regularization of the Laplacian determinant and the Fredholm de- 
terminant (|4.11|l are closely related. In Ref. (2Sj a general relation between them 
is found for any sequence of numbers of arbitrary order ^. In our particular case 
this relation reads 



dZ(s,0) 



ds 



dZ(s,a) 



s=0 



111 



s=0 



n 1 



,fc=0 



•oFP[Z(l)] (4.16) 



where FP[Z(1)] denotes the finite part of Z at 1, defined as usual by 



¥P[Z{ 8 )\ = 



'Z(s) 

FP[Z( S )} = Urn 

e— >0 



Z(s + e) 



ResZ(s) 



if s is not a pole of Z 
if s is a pole of Z 

(4.17) 



IV. 2 The Connection with the Grand Potential of the 
Coulomb Systems 

As the reader probably noticed, the Fredholm determinant is quite similar to 
the infinite product that appears in the expression for the grand canonical par- 
tition function (|2.5|l for Coulomb systems in the Debye-Hiickel approximation. 
Our goal in this section is to relate this kind of products with the geometrical 
information that can be extracted from the asymptotic expansion of the heat 
kernel. 



IV.2.1 The bulk case 

First let us mention some points concerning the case of an unconfined system. 
In this case the eigenvalues A n = A^ and the expression 1|2.5[) for the grand 
canonical partition function involves precisely the Fredholm determinant l|4.11|l 
with a = k 2 . The exponential terms e K / A ™ that come from the subtraction of 
the self-energy properly regularize the infinite product — (k 2 /A„)]. The 

final result for the grand potential l|2.6[l and (|2.7[) is finite and does not depend 
on any ultraviolet cutoff. 

As a side note let us mention that if we were interested in the formulation 
of Debye-Hiickel theory for a system living in four or more dimensions, the 
expression ()2.5(l for the bulk grand partition function would not be convergent 
and it would require the introduction of an ultraviolet cutoff. This is because 
in dimension d ^ 4 the index of the sequence of the Laplacian eigenvalues 
would be^ = d/2^2. For this case the correct regularization of the Fred- 
holm determinant would require an additional exponential term as shown in 
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Eq. l|4.1 (J [I . Physically this means that the bulk properties of a Coulomb system 
in dimension greater or equal than four, in the Debye-Hiickel regime, can only 
be defined for a system of charged hard spheres or any other charged particles 
with a short-range potential that cuts the singularity of the Coulomb potential. 
The inverse of the radius of the particles is the equivalent of the ultraviolet cut- 
off in our formulation. The bulk thermodynamic properties would depend on 
the radius of the particles, and diverge if one takes this radius to zero. This can 
be contrasted with the two- and three-dimensional cases where one can build a 
Debye-Hiickel theory for which the bulk properties have a well defined limit for 
point-like particles. 

As a complement on this remark let us remind the reader that for a three- 
dimensional system the exact thermodynamic properties, beyond the Debye- 
Hiickel approximation, are not well defined for a system of point-like particles 
due to the collapse of particles of opposite sign. In two dimensions this collapse 
problem is less strong: if the thermal energy of the particles is high enough a 
system of point particles is well defined. On the other hand, the Debye-Hiickel 
approximation is less sensitive to this collapse problem: for two- and three- 
dimensional systems the bulk properties are well defined for point particles. 
However as we have seen in the examples, in the three-dimensional case the 
surface properties are sensitive to the collapse problem and a proper definition 
of these require the introduction of a short-distance cutoff. For dimensions 
greater or equal than four the collapse problem appears for the bulk properties 
even in the Debye-Hiickel approximation. 

IV. 2. 2 Zeta regularized grand potential 

Now let us consider a confined Coulomb system in the manifold Ai g . Let R 
be the characteristic size of the manifold. For instance one can define R as 
y x l d where V is the volume of the manifold. We are interested in the large- R 
expansion of the grand potential of the system, which can be obtained from 
Eq. H2.5fl . In this section we will study an intermediary quantity £1* related to 
the grand potential which we define by 



For obvious reasons (remember Eqs. (|4. 1 i5|) and 1)2. 5|l ) we will call this quan- 
tity the zeta regularized grand potential. The prime in Eq. (|4.18|) indicates 
differentiation with respect to the first variable of the zeta function (s). 

As we will show below, the large- R expansion of the zeta regularized grand 
potential is related to the small-t expansion of the heat kernel (|4.3|) . To see this, 
let us consider a system where all lengths have been rescaled by a factor 1 /R : 
it is a Coulomb system confined in a manifold of fixed volume equal to 1 but 
with the same shape as the original system. Let Zi(s, a) be the generalized zeta 
function for the spectrum of the Laplace operator for such a manifold and 0i(i) 
its heat kernel. The subscript 1 refers to a system confined in a volume 1. For 




(4.18) 
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the original system of characteristic size R we will eventually use the subscript 
R in the spectral functions Zr(s,o) and &n(t). 

The eigenvalues of the system of size R are the same as those of the system 
of size 1 multiplied by a factor R~ 2 . Then we have Qn(t) = Q% (R 2 t) and 
Zn(s,a) = R 2s Zi(s,aR 2 ). From this we see that an expansion of the heat 
kernel Qr for large- R and fixed argument is the same as an expansion of the 
heat kernel 0i for small argument. 

We have 



§- s [R^ Zl (s,^R% 



(4.19) 



= Zi(0,0)lni?- Z 1 (0,K 2 R 2 )\nR- 



Z[(0,k 2 R 2 ) Z[(0,0) 



where the prime denotes differentiation with respect to the first argument s. The 
last term is a constant, so we will eventually drop it in the large- R expansion. 

From Eq. 14.1 2 bfl we can see that a small argument t expansion of the heat 
kernel Q(t) is equivalent to a large argument a expansion of the zeta function 
Z(s,a). Then using the small argument expansion of the heat kernel (|4.2() into 
Eq. d4~T2hll one can obtain in general the large- R expansion j2H] 



Z[(0,( K R) 2 ) 



T(i n )( K R)- 



m=0 



In(nR) 2 -Vr 1 



r=l 



(-K 2 R 2 y 



E * 

i„£z-u{o} 

(4.20) 

In this equation the coefficients Ci n are those of the heat kernel expansion (|4.2|l 
for a system of size 1 and we use the convention that if m is not any of the 
exponents i n of Eq. (|4.2|) then c rn = 0. The first sum runs over all indexes i n 
that are not negative or zero integers. Using this equation and [21 Zi(0,a) — 
c — ac_i into (|4. 19|> yield the large- R expansion of the zeta regularized grand 
potential 



R — > oc 



c ln T{i n ){KR)~ 



k 2 R 2 c- X 



In k 

2 



+ c \ii(kR) + -Z'(0,0) (4.21) 



Now we specialize this result for two and three dimensions. In two dimensions 
from Eq. (|4.9|) we read c_i = V/(4tt), c_ 1/2 = -B/(8^/n) and c = x/6 
where V = V/R 2 and B = B/R denote the volume (area) and perimeter of the 
manifold of characteristic size 1. Then we have 



m*2D = 



n 2 R 2 
Air 



1 _ l n K \y - + | \u(kR) + 0(1). (4.22) 



In three dimensions from Eq. I|4.3|l we have c_ 3 / 2 = V / (Att) 3 / 2 , c_i = 
-B/(16tt) and c_ 1/2 = (2(7 + j D)/[6(4tt) 3 / 2 ] with V the volume of the sys- 
tem of size l, B the area of the boundary of system of size 1 and C and D 
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the curvatura integra and integrated mean curvature for the system of size 1. 
Replacing this into Eq. (|4.21() we have 



• 2z?2 r 



pn* 3D = - 



K R ^ I , ' ft 



12tt 



16tt 



In k — 



B + ^(2C + D)+o(R) (4.23) 



IV. 2. 3 Connection between the physical grand potential and the 
zeta regularized grand potential 

The excess grand potential £l oxc of the Coulomb system is obtained from Eq (|2.5|l 

as 



1 



Mn 1 



K 
An 



(4.24) 



This expression involves a product very similar to the Fredholm determinant (|4.11|) 
but it only coincides with it for the bulk case. In general for a confined sys- 
tem they are different. However we can formally make appear the Fredholm 
determinant F by writing 



/30 e 



1 



HIT 



K 

Aw 



E, 



-E, 



2 kf M + T 



Ezrr-E 



-AS 



(4.25) 



Then we make the connection with the zeta regularized grand potential Q* 
defined in the previous section using Eq. I|4.16|) . so finally 



K 2 K 2 

JLfp[z(i)] + ^ 



_ \ X] _\0 



(4.26) 



The last two sums are divergent when taken separately. In principle they should 
be cutoff in a similar way as in the examples. The proper treatment of these 
sums should be done separately for each dimension d = 2 and d = 3. 



IV. 2. 4 Two-dimensional case 

In two dimensions the divergences of the sums involved in Eq. ()4.26(l can be 
dealt in an elegant way by means of the zeta function. The zeta function has 
a pole at s = 1. Remembering Eq. (|4.13(l and the fact that the residue of the 
zeta function Z at s = 1 for our confined system is equal to the residue of Z° 
at s = 1 for an unconfined system we can identify the summations in Eq. I|4.26|l 
with 

Am s— 
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Then using the fact that both zeta functions Z and Z° have the same residue 
at s = 1 we find that Eq. (|4.26() becomes 

pnfs = pn* D - ^ fp[z°(i)] . (4.28) 

Now the zeta function for an unconfined system reads 

z ° (s) = ? F¥ = v Lw?^ «« 9 » 

However this zeta function cannot be properly defined: if 5fte(s) ^ 1 the integral 
is ultraviolet divergent and if 5Re(s) > 1 it is infrared divergent. Depending on 
the sign of 5fte(s) — 1 this zeta function should be regularized with an ultraviolet 
or infrared cutoff. For our present purposes we need Z{s) defined for 5Re(s) > 1, 
then introducing the infrared cutoff k m i n = 2e~ c / L as in the two dimensional 
examples of section ITTT1 and Ref. JHj we have 

Z ( s ) = -H^sis- (4.30) 
w 4?r s - 1 

Its finite part at s = 1 is 

FP[Z°(1)] - -^lnfc min - -illn^. (4.31) 

Z7T Z7T L 

f3n?D = W*2n + ^ In ^ • (4.32) 

Clearly this extra term only contributes to the bulk grand potential. Finally 
reporting this into Eq. I|4.22|) yields the finite-size expansion 

^Td = ^ Q - C - In ^ ) V - + | ln( K R) + 0(1) . (4.33) 

We have obtained the general finite-size expansion of the grand potential for 
arbitrary confining geometry in two dimensions. We recover from the first term 
of Eq. I|4.33(l the bulk grand potential (|2.6() , from the second term the surface 
tension already obtained in the examples 7 = ksTn/8 and finally the logarith- 
mic finite-size correction (x/6) lni?. This constitutes a proof of the existence of 
this universal finite-size correction for Coulomb systems, in the Debye-Huckel 
regime, confined in an arbitrary geometry with Dirichlet boundary conditions 
for the electric potential. 

IV. 2. 5 Three-dimensional case 

For a three-dimensional system we must proceed differently to evaluate the sums 
in Eq. I|4.26[) as it was done in the two-dimensional case. Here we cannot identify 
the sums with zeta functions because in three dimensions the definition 14. 12a.ll 



Therefore 
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of the zeta functions expressed as a sum is only valid for 9?e(s) > 3/2 and in 
Eq. (|4.26() the sums correspond to s = 1. In the two-dimensional case we did 
not have this problem because the validity of l|4.12a|l if for 5ie(s) > 1 and since 
the residues of both zeta functions are equal we could take the limit s — > 1 + 
of the difference of zeta functions and obtain a finite result. Now, in the three- 
dimensional case, we need the sums for a value s = 1 < 3/2 which is far beyond 
the validity of Eq. (|4.12a() . Also the corresponding analytic continuations of the 
zeta functions do not have the same residue at s = 1. Indeed the bulk zeta 
function Z° does not have a pole at s — 1 in the three-dimensional case but the 
zeta function Z for the confined system has a pole at s = 1 with residue given 
by Eq. H4.13fl which is related to the coefficient c_i corresponding to the surface 
contribution to the grand potential. This suggest that the difference of the two 
sums in Eq. (|4.26|l will not be convergent for the three-dimensional case as it 
was in the two-dimensional one and it would give a surface contribution to the 
grand potential. 

Let us introduce a truncated version of the zeta function evaluated at s = 1 , 

z cut (x)= E A- ( 4 - 34 ) 

|A fc |<A 

and the corresponding one for the unconfincd system 

I A fc |<A 

Here A > is an ultraviolet cutoff for the eigenvalues. The sums in Eq. 14.26fl 

are 

E A~ - E Ao = - lim [^ut(A) - Z cut (\)] . (4.36) 

m n b 

Let us introduce the counting function which is equal to the number of 

eigenvalues Afe which are |Afc| < A. The truncated zeta function is related to the 
counting function by 

Zcut(A)= t^AdX (4.37) 
Jo A 

with Af'(X) the derivative (in the sense of the distributions) of jV(A). On the 
other hand the derivative of the counting function and the heat kernel are related 
by a Laplace transform 

/>oo 

0(t) = / e - tx Af'(\)d\ (4.38) 



Then, from the asymptotic expansion of the heat kernel l|4.2(l for t — > 0, we can 
obtain the first terms of the asymptotic expansion of AT'(\) for A — > +oo 

JV(A)_~ Twfe^ 1/2 + c - 1 + ( 1 ) = 7^^ 1/2 -^- + °( 1 ) ( 439 ) 
A^oo 1 (3/2) (2ny lbn 
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Reporting this into Eq. (|4.37|) we find the asymptotic behavior of the truncated 
zeta function for A — > oo 

^cut(A) = -^A 1 / 2 --|-lnA + 0(l) (4.40) 

Then we have 

Z cu t(A) - Z c ° ut (A) = ~lnA+ 0(1) = In A max + 0(1) (4.41) 

107T 07T 

where we introduced A max — A 1 / 2 to make the connection with the examples 
of section HTT1 We see that the divergences are not completely canceled now as 
opposed to the two-dimensional case. There remains an ultraviolet divergence 
which contributes to the surface tension (proportional to B the area of the 
boundary of the manifold). This is the same kind of divergence that we have 
found in the three-dimensional examples in section lTTTl and it is due to the strong 
attraction of the particles to their images as explained earlier. 

Finally collecting the results from Eqs. (|4.23() . (|4.26() . and 14.41[> the finite- 
size expansion of the grand potential reads 



K 3n3 2d2 

A 3D 12tt 16tt 



In-— + 0((A max ) ) 



B+—(2C + D)+o(R). 

487T 

(4.42) 

From this very general calculation we recover the bulk grand potential Il2.7|) . 
the surface tension 7 = fcsT(«; 2 /167r) ln[/t/if max ] + 0(1) which is ultraviolet 
divergent. Also we find the perimeter correction to the grand potential (the 
term proportional to R) which depends on the curvatura integra C and the 
integrated mean curvature D. 

In the example of section Unl of the three-dimensional Coulomb system con- 
fined in a ball we have found a perimeter correction to (3Q equal to kR/3. 
For a ball of radius one in the flat space R 3 the curvatura integra is C — 
and the mean curvature computed from Eq. I|4.8|l is J = 4 and the integrated 
mean curvature is D — 16ir. Then the correction predicted by Eq. (|4.42() . 
kR(2C + £>)/(487r) = kR/3, is in agreement with the explicit result found in 
the example. For the other example of the Coulomb system confined in the 
thick spherical shell, the agreement of our general result (|4.42(1 with the explicit 
calculation of section I1TTI is also straightforward. 



IV. 3 System Confined in a Square Domain 

The above analysis is valid for a confining manifold with smooth boundary. How- 
ever it can easily be generalized for a manifold with a boundary with corners. 
As an illustration let us consider the case of a two-dimensional Coulomb system 
confined in an square domain of side R and subjected to Dirichlct boundary 
conditions for the electric potential. 

Conformal field theory predicts that in the case of a two-dimensional crit- 
ical system confined in a geometry with corners in the boundary, it appears 



22 



a contribution to the free energy (times the inverse temperature) equal to 
[6/ (247r)](l — (tt/9) 2 ) In R for each corner with interior angle 6 [7j- In the case of 
a square 9 = n/2 and the contribution per corner equals fjf (1 — (Stt/tt) 2 ) lni? = 
— (1/16) lni?. Then the total contribution of the four corners is —(1/4) lni?. If 
the similarity of conducting Coulomb systems with critical systems holds then 
we should expect a finite-size correction to the grand potential times (3 equal to 
(1/4) lni? for Coulomb systems. 

The eigenvalues for this case can be found easily by separation of variables 
expressing the Laplace operator in rectangular coordinates and solving the eigen- 
value equation. The spectrum for a system in a square of side equal to 1 is given 



byA n>i = -7r 2 (n 2 + Z 2 ), 



1,2,... and 1= 1,2, . 



The heat kernel is 



k=0 



which can be expressed in terms of the Jacobi theta function 



(4.43) 



(4.44) 
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as 



e 1 (t) = 1 [0 3 (om) - if 



(4.45) 



The heat kernel expansion for t — > can be found using Jacobi imaginary 
transformation |31) 



This gives 



ei(t) = - A 



-l + iirt)- 1 / 2 l + 2^e-" ; 



/* 



(4.46) 



(4.47) 



The asymptotic expansion for t — > is 



1 



i i 



4?rt 2\/7rf 4 



1 + 



,-i/t 
t 



(4- 



Comparing with expression l|4.9|) for a smooth boundary we recognize the first 
two terms: the volume (area) (V — 1) and the surface (perimeter) (B = 4) 
terms which are the same. The constant term on the other hand is now equal to 
1/4. Applying the same argument developed above for the general case to this 
heat kernel we see that this constant term is the one that gives the coefficient of 
the logarithmic finite-size correction for the grand potential. Then the finite-size 
expansion for this geometry reads 



/?fi = (3tt b - — + - In(nR) + o(ln R) 
8 4 



(4.49) 
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with Hi, given by Eq. (|2.ti|) with V = R 2 and the perimeter of the square B = 4R 
with R the length of a side. The logarithmic finite-size correction is in agreement 
with the one predicted by conformal field theory with the appropriate change 
of sign. 

V Summary and Perspectives 

We have illustrated with several examples how to apply the method of Ref. ^Bj 
to find the grand potential of a Coulomb system in the low coupling regime and 
confined by ideal conductor boundaries. We considered several examples: in 
two dimensions the slab and in three dimensions the slab, the ball and a thick 
spherical shell. The method can easily be adapted to other geometries. In all 
the examples we also computed the finite-size expansion of the grand potential. 
For the slab geometries, in three and two dimensions, we recover a universal 
algebraic finite-size correction predicted in Ref. P^. For two-dimensional fully 
confined systems the finite-size expansion exhibits a universal logarithmic term 
similar to the one predicted for critical systems by conformal field theory |16j. 

We have also extended the method to confined systems of arbitrary shape 
with a smooth boundary. For this general case we showed how the heat kernel 
expansion for small argument for the considered geometry is related to the large- 
size expansion of the grand potential of the Coulomb system. From this, we 
recovered the expressions for the bulk grand potential and the surface tension 
which agree with those found in the specific examples. Regarding the finite- 
size corrections, in the case of two dimensions we proved the existence of a 
universal logarithmic finite-size correction for the grand potential times (3 equal 
to (x/6)lni? with x the Euler characteristic of the confining manifold. For 
three dimensional systems we also found a general prediction for the perimeter 
correction to the grand potential but it is not universal (it depends on the Debye 
length) . 

The general treatment for arbitrary confining geometry exposed here is done 
for Dirichlet boundary conditions but can also be adapted for other kind of 
boundary conditions. For example for ideal dielectric boundary conditions, 
i.e. Neumann boundary conditions for the electric potential, the heat kernel 
expansion is very similar to Eq. (|4.3|l except that the surface term changes of 
sign ^21- As a consequence the surface tension of a Coulomb system confined 
by ideal dielectric boundaries would be minus the surface tension of a Coulomb 
system confined by ideal conductor boundaries in the low coupling regime. For 
a confining geometry with smooth boundary the finite-size expansion will read 
as in Eqs. I|4.33[) and 14.42(1 for two and three dimensions respectively except 
for a change of sign in the surface tension term. 

The method exposed in Ref. ^Hj an d used here can be extended to compute 
the density profiles and correlation functions. An example of such application 
for the case of a slab geometry in three dimensions can be found in Ref. |19) . 
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